'Algorithmic cooling' in a momentum state quantum computer 
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We describe a quantum computer based upon the coherent manipulation of two-level atoms be- 
tween discrete one-dimensional momentum states. Combinations of short laser pulses with ki- 
netic energy dependent free phase evolution can perform the logical invert, exchange, CNOT and 
Hadamard operations on any qubits in the binary representation of the momentum state, as well 
as conditional phase inversion. These allow a binary right-rotation, which halves the momentum 
distribution in a single coherent process. Fields for the coherent control of atomic momenta may 
thus be designed as quantum algorithms. 

PACS numbers: 03.67.Lx, 32.80.Pj, 33.80.Ps, 39.20.+q 



Proposed schemes for quantum computation [jjj have 
tended, quite naturally, to focus on quantum analogues of 
classical binary computing elements. The nuclear spins 
of a molecule, or of an ensemble of trapped atoms or 
ions, thus mimic the bits of a conventional computer. 
In this article, we address a less obvious system, in 
which information is represented by the momentum of 
a single atom or molecule, which is manipulated using 
laser pulses in a one-dimensional geometry that restricts 
each species to a ladder of equally spaced momentum 
states. Although any one laser pulse can change the 
species momentum, through photon absorption and stim- 
ulated emission, by only a single photon impulse, we find 
that sequences of pulses, interspersed with periods of 
momentum-dependent phase evolution, allow a full suite 
of quantum computational operations on the qubits com- 
prising the binary representation of the momentum state. 

The size of the momentum state quantum computer 
grows in proportion to the number of quantum states in- 
cluded, where conventional candidates instead scale with 
the number of qubits representing those states. The 
number of laser pulses needed to perform each logical 
operation increases similarly and, although the overall 
duration proves to be less drastically affected, momen- 
tum state systems therefore hold limited promise for real 
computing. Nonetheless, the scheme outlined here is 
based upon simple and readily-available elements, albeit 
in complex combinations, and could thus complement 
NMR systems || || as a testbed for experimental studies. 

It is, however, in the design of complex fields for co- 
herent control that we forsee the greatest potential, for 
if momentum-changing operations can form the basis of 
a quantum computer then the pulse sequences for op- 
tical manipulation may be optimized as quantum com- 
putational algorithms. In this respect, the momentum 
state quantum computer is an enthusiastic development 
of schemes for interferometric cooling B and the coher- 
ent amplification of laser cooling f| . 



Our scheme is based upon the motion in one dimension 
(henceforth taken to be vertical) of a sample of two-level 
atoms, such as an atomic beam interacting with trans- 
versely propagating laser beams, as shown in Fig. |l|. Wc 
shall refer to four coherent operations: 

W+(a,4>) a short upward laser pulse 
VF_ (a, 4>) a short downward laser pulse 
F{ijjt) free evolution (electronic energy) 
G(t/r) free evolution (kinetic energy), 

where r = 2m/(hk 2 ). Here, the short (and therefore 
spectrally broad) laser pulses couple the upper and lower 
atomic levels, between which population is transferred 
through Rabi cycling for the duration of the pulse. Con- 
ventionally, we describe the overall effect of the pulse 
through the phase 2a of the Rabi cycle incurred, the 
population being inverted when 2a = it (the so-called 
f 7T pulse'), and restored when 2a — 2n. Other fractions 
of a 7r pulse will convert an initially pure state into a 
superposition. The phases <j> are determined by the rela- 
tive optical phases of the laser pulses. The free evolution 
operations F(uit) and G(t/r) correspond simply to the 
components of the time-dependent wavefunction phase 
exp(— iEt/h) that correspond to the electronic energy 
and vertical momentum component respectively. Weitz 
and Hansch pj have shown how the electronic and kinetic 
energy contributions to the free phase evolution may be 
separated by appropriate insertion of pairs of 7r-pulses 
that invert the atomic population, as we discuss later. 

In this one-dimensional geometry, the atom is con- 
strained to a ladder of momentum states that are spaced 
at intervals of the photon momentum Tik = Two j c (w be- 
ing the frequency of the resonant transition) and alter- 
nate between the ground and excited electronic levels g 
and e. We label these states according to their momen- 
tum components, in units of the photon impulse hk. We 
initially assume, for our analysis of the momentum state 
quantum computer, that these momenta take integer val- 
ues, but this assumption will be relaxed when we later 
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FIG. 1: (a) Fractional ir pulses tuned to the two-level atom 
couple adjacent momentum states, which (b) we label in units 
of hk. This number, in binary representation, gives the qubits 
of the momentum state quantum computer, (c) Two n/2 
pulses may act as the beamsplitters of an atomic interfer- 
ometer; the relative phase between the two paths determines 
whether the pulses add or subtract, and hence whether or not 
the electronic state is inverted. 



consider the consequences for atomic manipulation. 

We now convert our representation of the ladder of 
states to binary, using the notation Q n . . . QzQiQq, with 
the least significant bit on the right. This is the crucial 
step in our analysis; yet, apart from the correspondence 
of Qo to the electronic state of the atom, binary rep- 
resentation seems at first rather unpromising, for com- 
putational notation usually helps only when the bits 
themselves can be manipulated. While the momentum- 
changing laser pulses here move population by at most 
one state at a time, however, appropriate combinations 
prove to offer exactly such bit-wise manipulation. 

The key, as in the interferometric cooling scheme of 
Ref. jl} , is the dependence of the phase of free evolution 
upon the momentum. For two levels p + Ap/2 and p — 
Ap/2 (in units of hk) and electronic energy difference 
E 2 i = huj, the relative phase ip evolves according to 
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FIG. 2: Bloch-vector representation of the first stage (Eq. (g)) 
of the right-rotation. The first pulse rotates the four ground 
states into the horizontal plane; free evolution distributes 
these around the vertical axis according to their momenta; 
the second pulse then returns two to pure states, leaving the 
others in mixed states. Aside from phase corrections, the full 
right-rotation operation takes 18 n/2 pulses and 26 n pulses. 



cause inversion (the n pulse) and when 2a = 0, 2ir and 

so on. For periods of free evolution, we simply indicate 
the relative phases introduced between coupled states. 

First consider a pair of upward-travelling tt/2 pulses, 
separated by a period of free evolution that introduces 
between coupled states a relative phase adjusted to give 
in each case an integer power of exp(— i7r/2). This forms a 
simple interferometer which, depending upon the original 
state momentum, can invert population, return it to its 
original state, or leave initially pure states coupled (see 
Fig. 0(c)). The sequence repeats every 8hk. 
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Any pair of momentum states thus incur a relative phase 
that evolves according to their average momentum p. 
Cancellation of the electronic contribution to the phase, 
by inserting n pulse pairs so that the states spend equal 
times in the ground and excited levels [Q , merely changes 
the average Ap and hence the rate at which this proceeds. 

We illustrate the capacity for bit-wise manipulation 
with the example of a three qubit right rotation, 

{Q2,Qi,Qo} —> {Qo,Q2,Qi}- 

In our largely diagrammatic description, which has its 
origins in Fig. [l](c), pulses or complete pulse sequences 
coupling adjacent levels are indicated by ^ when they 
produce a superposition (e.g., a 7r/2 pulse), X when they 



This and subsequent operations may be visualized, as 
in Fig. |2|, as rotations of Bloch vectors representing the 
coupled states - a picture also suitable for NMR com- 
puters ||, Pure states are vertically up (\g)) or down 
(|e)), and are coupled by rotation about a horizontal axis, 
whose direction depends upon the optical phase. Free 
evolution corresponds to rotation about the vertical axis. 

Two such sequences may now be combined with a fur- 
ther momentum-dependent phase between them to form 
a second interferometer for the states left in superposi- 
tions. The ladder of phases is offset either via the F{uit) 
operation (such as a period of uncompensated free evolu- 
tion that is dominated by the 'electronic' phase ), or by 
appropriate phasing of the subsequent laser pulses. The 



result is a conditional state exchange: 
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Precisely which pair of adjacent states is exchanged by 
this operation depends upon the directions and relative 
phases of the tt/2 pulses. As the penultimate step, we 
construct the two-qubit exchange operation EX(2,1): 
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The right-rotation RR3 is completed by combining this 
sequence with the much simpler exchange EX(1,0) of 
qubits Qi and Q . This is given in Table | together 
with a range of other operations on the first three qubits 
from which, in conjunction with the one-bit operations 
W + (a, 4>) and F(ojt), a complete set may be formed ||. 

When only the ground electronic level is occupied 
(Qo = 0), the right-rotation is indistinguishable from a 
divide-by-two operation (e.g., 100(= 4) — > 010(= 2)) and 
provides a cooling mechanism. As shown in Fig. |3| the 
four ground states are coherently mapped onto the four 
lowest momentum states, two of which subsequently un- 
dergo spontaneous emission, leaving only states 0, 2 and 
4 populated. The coherent process, which pumps heat 
from kinetic to electronic energy, may then be repeated, 
further narrowing the momentum distribution. 

At this point, we relax our earlier assumption of 
integer-valued momentum states, and find that the ef- 
fect, whilst imperfect for non- integer momenta, nonethe- 
less remains. The results of our simulations for the effect 
of this sequence on an initially flat distribution across 
fractional momentum states are shown in Fig. ^. 

For our simulations, we have used matrix representa- 
tions of the pulse and evolution operations. Although the 
matrices are in principle infinite, all non-zero terms clus- 
ter around the leading diagonal and any element rriij dif- 
fers from the diagonally displaced term mi + 2 n ,j+2n only 
through its momentum dependence, so we may summa- 
rize the matrices as 4 x 4 elements, given below and de- 
rived from the equations of Friedberg and Hartmann . 
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FIG. 3: Cooling via the right-rotation operation, shown here 
applied to the lowest three qubits. The initial distribution 
across ground states {0, 2, 4, 6} is transferred (bold arrows) 
to the lowest momentum states {0, 1, 2, 3}; subsequent spon- 
taneous emission leaves population in states {0,2,4}. The 
width of the momentum distribution may thus be reduced by 
a factor of nearly 2 in a single coherent step. 
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FIG. 4: Simulated evolution of the momentum state distri- 
bution, shown after 1, 2, 4 and 8 cycles of the 3-qubit coher- 
ent cooling algorithm. Although the process applies perfectly 
only to atoms with even, integral momenta, significant cooling 
remains apparent. Spontaneous emission scrambles the exact 
momenta. After only a few coherent cycles, the distribution 
has been narrowed to less than a single photon impulse. The 
initial probability density is unity. 



Matrices that we use in practice must merely be ex- 
panded to cover the sequence of interactions and the mo- 
mentum range that we wish to describe. The following 
matrices act on the states {2, 1, 0, —1}. 

For upward and downward travelling fractional 7r 
pulses corresponding to Bloch vector rotation through 
'Rabi angle' 2a and optical phase </>, we have 



W+(a,4>) = 



( cos a \ 
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(5) 



and 
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ie^sina cosa 
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y ie 1 ^ sin a cos a J 



(6) 

The matrices for free evolution according to the electronic 
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the large impulse per coherent cycle would be a particular 
advantage. 



(7) 



and 



( e -i(Po+2) 2 t/r 

e -i( P0 + l) 2 t/T Q 

e-iPot/r o 

V e -i(po-i) 2 t/r J 

(8) 

where r = 2m/(hk 2 ). Weitz and Hansen's sequence for 
interferometric cooling thus becomes 



F(2cT) • ' G(^) ■ F{ljT') ■ W + Q 



No attempt has yet been made to optimize the se- 
quences given here: the phase corrections often serve 
only to demonstrate the exact equivalence to a quan- 
tum computer, and reductions in the complexity, dura- 
tion or momentum sensitivity of each operation should 
be possible. Nor have we examined superpositions in- 
volving more than two states flio| or interactions at more 
than one wavelength . Instead of the fractional Rabi 
coupling of electronic transitions assumed for simplic- 
ity, our scheme could be more robustly implemented us- 
ing Raman transitions for adiabatic passage (l^, Q 
between Zeeman or hypcrfmc levels, possible even with 
modulated c.w. lasers The scheme, which we think 
of as a form of 'algorithmic cooling' |l6| in its broadest 
sense, could in principle be extended to three dimensions. 
Owing to the non-resonant nature of the pulsed interac- 
tions, it would also be suitable for molecules, for which 
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TABLE I: Basic operations of the momentum state quantum computer. No attempt has been made to optimize the pulse 
sequences, which run from right to left, and some uncorrected phases remain in the operations marked with an asterisk. For 
the operation G, po is taken to be zero (mod 8) . 



level 


name 


description 


sequence 


basic 


G(t/r) 




W- (tt, 0) • FG(t /4t) ■ W- (tt, 0) • FG(t/4r) ■ 
W+(ty, 0) • FG(t/4r) ■ W+(n, 0) • FG{t/Ar) 


1 qubit 


NOT(O) 


Qo->Qo 


^(f)-^+(!>o)-F(f) 




CP1(0) 


if state=0, invert phase 


F(n)-W+(n,0) 




HAD(O) 


Walsh-Hadamard on Qo 


W+($,%)-F{n)-W+(*,0) 


2 qubit 


EX(1,0) 


{Qu Qo} -» {Qo, Qi} 


^(f)-^-(f^)-G(f)-W_(f,f).F(^) 




CNOT(1,0) 


{Qi,Qo} -> {Qi,Qi0Qo} 


F(f).W/ + (f, 7 r).G(f).W + (f,f). J F(5f) 




CNOT(l.O) 


{QuQo} {Qi,Qi®Qo} 


M/+(^,0)-F(^) • W+(f,0) -G(f) • W+(f,f) -F(^) 




CP2(0) 


if state=0, invert phase 


F(^)-G(f).M/ + (7r,0) 




HAD(1,0) 


Walsh-Hadamard on Qi, Qo 


EX(1,0) • HAD(O) • EX(1,0) • HAD(O) 


3 qubit 


SW3(2,3)* 


swap states 2, 3 


W + (f,0)-G(f).W + (f,f)-F(f).G(f). 

w+(i>°) • G (i) • w +{h f ) • fffi • G(f ) 




SW3(3,4)* 


swap states 3, 4 


F(tt) • M/_ (f,0) • G(|) • , f ) • F(^) • G(f ) ■ 
W-(f,0) -G(f) -W-(f,f ) -G(f) 




SW3(4,5)* 


swap states 4, 5 


^+(l 1 0).G(f).^ + (f,f).F(^).G(f). 




EX(2,1) 


{Qa,Qi}-{Qi,Q a } 


H/ + (f,0).G(f).H/ + (f,f).F(f).G(f) 
H/+(tt,0) • CNOT(1,0) • EX(1,0) • G(^) ■ F(ift) ■ EX(1,0) ■ CNOT(1,0)- 
SW3(3,4) ■ NOT(0) ■ F(ir) ■ NOT(0) ■ SW3(4,5)- 
NOT(0) ■ F(tt) ■ NOT(0) ■ SW3(2,3) ■ SW3(3,4) ■ G(^f) ■ F(^) 




RR3 


{Qa,Qi,Qo}-{<?o,Qa,Qi} 


EX(2,1) • EX(1,0) 




RL3 


{Q 2 ,Qi,Q } {Qi,Qo,Q 2 } 


RR3 • RR3 




CP3(0) 


if state=0, invert phase 


NOT(0) • RL3 ■ NOT(0) ■ RL3 ■ F(^) • G(^) ■ RR3 ■ SW3(4, 5) ■ F(^) ■ 
SW3(4, 5) ■ F(f ) • NOT(0) • RR3 ■ NOT(0) ■ W+(n, 0) 



